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Abstract 

The Nambu~Jona-Lasinio (NJL) model has been extensively studied by many researchers. In 
previous work we have generalized the NJL model to include a covariant model of confinement. 
In the present work we consider further modification of the model so as to reproduce the type of 
Euclidean-space momentum-dependent quark mass values obtained in lattice simulations of QCD. 
This may be done by introducing a nonlocal interaction, while preserving the chiral symmetry of the 
Lagrangian. In other work on nonlocal models, by other researchers, the momentum dependence of 
the quark self-energy is directly related to the regularization scheme. In contrast, in our work, the 
regularization is independent of the nonlocality we introduce. It is of interest to note that the value 
of the condensate ratio, {ss) / {ud), is about 1.7 when evaluated using chiral perturbation theory and 
is only about 1.1 in standard applications of the NJL model. We find that our nonlocal model can 
reproduce the larger value of the condensate ratio when reasonable values are used for the strength 
of the 't Hooft interaction. (In an earlier study of the r;(547) and r/'(958) mesons, we found that 
use of the larger value of the condensate ratio led to a very good fit to the mixing angles and decay 
constants of these mesons.) We also study the density dependence of both the quark condensate 
and the momentum-dependent quark mass values. Without the addition of new parameters, we 
reproduce the density dependence of the condensate given by a well-known model-independent 
expression valid for small baryon density. The generalization of our model to include a model of 
confinement required the introduction of an additional parameter. The further generalization to 
obtain a nonlocal model also requires additional parameters. However, we believe our results are 
of sufficient interest so as to compensate for the introduction of the additional parameters in our 
formalism. 
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I. INTRODUCTION 



The Nambu-Jona-Lasinio model has been extensively studied for several decades [1-3]. 
In recent years there has been strong interest in the study of quark matter at high densities, 
using the NJL model and related models. In particular, one finds color superconductivity 
under certain conditions and it has been suggested that some compact stars might be made 
of superconducting quark matter [4-10]. It is our belief that in such studies one should use 
a model which reproduces, as well as possible, know features of QCD. In this work we wish 
to generalize the SU(3)-flavor version of the NJL model to be consistent with the type of 
momentum-dependent quark masses found in lattice simulations of QCD [11]. For example, 
in Figs. 1 and 2 we show some of the results obtained in Ref. [11]. It may be seen from 
these figures that, in Euclidean space, the quark mass goes over to the current mass for 
Euclidean momentum k > 2 GeV. On the other hand, the NJL model, in the standard 
analysis [1], gives rise to a constant value for the constituent quark mass. As we will see, 
the form of the nonlocahty used here is different from that used in Refs. [6, 12-16]. For 
example, in reference [12] the qq vertex is modified by a form factor which depends on the 
relative momentum of the quark and antiquark, while, in another scheme, a form factor is 
associated with each quark line appearing in a diagram. In the latter procedure no further 
regularization is needed. However, for the problem considered in this work, the nonlocal 
models that appear in the literature are of limited applicability, since, in the limit of zero 
current quark mass, the quark self-energy is proportional to the form factors used to define 
the nonlocality [12]. In contrast, in the current work, we calculate the form of the quark 
self-energy after introducing a regulator and a nonlocal quark interaction. We stress that 
the procedure used here differs from any that appears in the literature. 

The organization of our work is as follows. In Section II we review the standard analysis 
for the condensates and "gap equation" of the SU(3)-fiavor NJL model [1]. We then go 
on to review a procedure for including the contribution to the quark self-energy due to the 
addition of a confining interaction for the case of the SU(2)-flavor model. In Section III 
we introduce a nonlocal interaction in the SU(3)-flavor NJL model while maintaining the 
separable nature of the interaction. We also describe the approximation used for the 't Hooft 
interaction in the nonlocal model. In Section IV and V we present some of the results of 
our numerical calculations. In Section VI we discuss the density dependence of the quark 
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FIG. 1: Quark mass values obtained in Ref. [11] for various current quark masses: mP = 91 MeV 
[ circles] , mP = 54 MeV [crosses] and mP = 35 MeV [diamonds] . 



condensate and the momentum dependent quark mass. Finally, Section VII contains some 
additional discussion and conclusions. 



II. THE QUARK SELF-ENERGY IN THE NJL MODEL 

In order to best introduce the nonlocal model, we will first review the calculation of the 
quark self-energy in the local SU(3)-flavor NJL model and then proceed to add a confinement 
interaction, as was done in an earlier study of the quark self-energy [17]. We consider the 
generalization to a model with nonlocal short-range and 't Hooft interactions in the next 
section. The Lagrangian of the model is 

jC = q(i^ - m')q + ^ Yl^iq^qf + {qh.yqf] 

+ ^{det[g(l + 75)?] + det[g(l - 75)?]} . (2.1) 

Here mP is the matrix of quark current masses, mP = diag (m°, m^, m°) and the A*(i = 
1, • • • ,8) are the Gell-Mann matrices. Further, Aq = •\/2/3 1, with Ibeing the unit matrix 
in the flavor space. 
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FIG. 2: Quark mass values obtained in Ref. [11] using an extrapolation of the current quark mass 
to zero. (The small dip at A; ~ 1.6 GeV is not statistically significant [11]. 



The quark propagator is written as 

iS{k) 

with 



We may define 



1 - Bu{k^) ' 



and 



(2.2) 



(2.3) 



(2.4) 



(2.5) 



1 - B^{k^) ' 

with similar definitions for Md{k'^), Ms{k'^), Zd{k'^) and Zs{k^). 

In the absence of a confinement model, we have B{k'^) — Q,Ay_{k'^) — Ad{k'^) — rUu and 
Ag{k'^) — nis, where m„ and and are constants. (Here, we have take = m^.) In this 



case we have [2] 

m„ = m° - 2Gs{uu) - GD{dd) (ss) , (2.6) 
md^m^a- 2Gs{dd) - Gd{uu) (ss) , (2.7) 

m, = m° - 2Gs{ss) - Gd{uu) (dd) . (2.8) 

These equations are depicted in Fig. 3a, where the last term represents the 't Hooft inter- 
action. The up quark vacuum condensate is given by 

d% C(P) ^ ^ 

= -4Nj^J"^f\ . (2.10) 

J (27r)4 k^-ml + ie ^ ^ 

Here, C{k'^) is a function needed to regulate the integral. In this work we will use the Pauli- 
Villars procedure and evaluate the integral in Euchdean space, as was done in Ref. [17]. In 

the general case we may write 

{UU) — —^ISA I ,^ —TT, , ..o/, ox -■ 5 (2.11) 

where Zu{k^) and M^{k'^) were defined in Eqs. (2.4) and (2.5). 

In the appendix of Ref. [17] we considered Lorentz- vector confinement, with 

V~-{kE - k'E) = r{'^)l^.{'^)V%kE - k'^) (2.12) 

and 

V%kE - k^) = -87r« |— ] - — ] , (2.13) 

where k^—k'^ denotes the Euclidean-space momentum transfer. Here, // is a small parameter 
introduced to soften the momentum-space singularities. Note that the form 



[ 1 V 



VHk - k') = -Sttac { ^ ) , (2.14) 

\[{k - k'Y + [{k - k'Y + ^ 



represents the Fourier transform of V'^{r) — Kre"'^'^, so that for small /i, y^(r) approximates 
a linear potential over the relevant range of r. In Fig. 3b we show the equation for the 
self-energy when the confining field is included. 
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In Ref. [17] we obtained the following coupled equations in the case of the SU(2)-fiavor 
model 

.n.^, _ . / [-AV^jk - k') + AN,njGs]A{k'^) 



These equations were solved after passing to Euclidean space and including a Pauli-Villars 
regulator of the form 

^^^^^ ^ [kl + A\kl) + K^][kl + A^{kl) + 2A2] ^^-^^^ 
in Euchdean space. Note that the form 

~ 2 2A^ 

^^^^^ ^ [kl (1 - + A^kl) + A2][^| (1 - B{kl)f + ^2(^2) + 2A2] ^^-^^^ 

may also be used. (In Eqs. (2.15) and (2.16) appears with a sign opposite to that given 
in Ref. [17], since, in that work, we used a negative value of k. For the present work we use 
a positive value of k to be consistent with all of our other publications.) 



III. A NONLOCAL NJL MODEL 

In this section we describe the procedures we use to create a nonlocal version of our 
generahzed NJL model. We consider the second term of Eq. (2.1) and make the replacement 

^ jyq{x)\,q{x)f + {q{x)i^,\,q{x)f] (3.1) 

i=0 

^ Y.{[q{x)yf{x)q{x) ■ q{y)X'mq{y) 

i=0 

+ [q{x)i'y5X'f{x)q{x) ■ q{y)il^\' f {y)q{y)\} ■ 

This replacement corresponds to the use of a separable interaction V{x — y) = Gsf{x)f{y). 

A related modification may be made for the 't Hooft interaction. It is useful, however, 
to describe these modifications as they affect momentum-space calculations. With reference 
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FIG. 3: a) A diagrammatic representation of the equation for the self-energy for a quark of momen- 
tum k. The first term on the right is the contribution of the current quark mass mP. The second 
term corresponds to the term proportional to Gs in Eqs. (3.6)-(3.8) and the last term represents 
the 't Hooft interaction. b)The self-energy equation in the presence of a confining interaction (wavy 
line.) Without the 't Hooft interaction, we have a representation of SU(2)-flavor model studied in 
Ref. [17]. (See Eqs. (2.15) and (2.16).) 

to Fig. 4, we replace Gs by f{ki — k2)Gsf{ks — ki). In the evaluation of the second term of 
Fig. 2 we need Gs{k — k') — f{k — k')Gsf{k — k') and choose to write 

Gs{k-k') = exp[-(A;-A;')'V2/3]G'5exp[-(A;-A;')'72/3] (3.2) 

= GseM-ik-k'Y"^/^]. 

In this work we take n — A and /3 20 GeV^. In Fig. 5 we exhibit the function Fik"^) — 
exp[— It is clear that many other functions may be chosen. 

Wc now rewrite Eqs. (2.15) and (2.16) for the up, down and strange quarks. For example, 
for the SU(3)-flavor case 

k^B (P\ - 2{k-k')[l-B^{k'')]V%k-k') 

^ ~ J (27r)4 fc'2[l - S„(fc'2)]2 - Al{k'^) + 16 ' ^ ^ ^ 

with similar equations for A^lk"^), Bdik"^), etc. Again, these equations are solved after passing 
to Euchdean space and introducing regulator functions: Cu{k'^), Cd{k'^) and Cs{k'^). In 
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iGs iGs(ki-k2,k3-k4) 



FIG. 4: The figure indicates the replacement of the local quark interaction, iGs, by the nonlocal 
(separable) term, iGs{ki — k2, ks — k4) = iGsf{k\ — k2)f{kz — k^C)- The distinction between our 
separable model and that of Ref. [12], for example, is that in Ref. [12] the alternative replacement 
iGs — > iGs{k2 + k4, ki + k^,) = iGsf{k2 + k4)f{ki + k^) was used. 




k (GeV) 



FIG. 5: The correlation function F{k) = exp[-A;^"/^] is shown for n = 4 and (3 = 20 GeV^. 
Euclidean space we write 

^"^^ ^ " + (P) + A2] [fc2 + (P) + 2A2] ' ^^"^^ 

etc. Note that without the 't Hooft interaction the equations for the up, down and strange 
quarks are uncoupled. 

Our treatment of the 't Hooft interaction is based upon a generalization of the last term 
in Eqs. (2.6)-(2.8). With reference to the third term on the right in Fig. 6, we introduce 
a correlation between the quark of momentum k and the quark of momentum k' . We 
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FIG. 6: The quark self-energy equation is depicted, with nonlocal terms replacing Gs and Gd- 
[See Fig. 4.] 



also include a correlation between the quark of momentum k and that of momentum k". 
(Therefore, our procedure does not introduce a correlation between the two quarks in the 
separate condensates. At this stage of the development of our model that seems to be a 
reasonable approximation and avoids having to define a three-quark correlation function, 

f{k — k',k' — k", k — k").) For example, we generalize the term —GD{dd){ss) to obtain the 
contribution to Au{k): 

dV Ca{k'^)Z,{k'^)M^{k'^)p{k-k') 



Kik) = -Gn 



-4:NJ 



X 



(27r)4 k'''[l - Bu{k'^)]^ - Al{k'^) 
dV Cs{k"^)Zs{k"^)Ms{k"^)P{k - k") 



(3.6) 



k"^[l- Bs{k"'^)Y - Al{k"^) 
in Minkowski space. This expression is then evaluated in Euclidean space and added to the 
right-hand side of Eq. (3.3). In a similar fashion, we calculate A^(^) ^s(^)- 



IV. NUMERICAL RESULTS: CONDENSATES AND CONSTITUENT MASS 
VALUES 

There is a good deal of flexibility in choosing the regulators C„(A;^), Cd{k'^)^ and Cs(/c^). 
Also, various forms could be chosen for the correlation functions, f{k — k'). Previously, in 
our Minkowski-space studies of the rj mesons we used Gs — 11.84 GeV~^ and Go — —200 
GeY~^ [18]. However, in that work we used a Gaussian regulator in Minkowski space so that 
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a direct comparison with the present study can not be made. On the other hand, we do not 
expect to find radically different parameters, if the constituent masses in the two calculations 
are similar. For example, in our earlier work, in which m„ and were parameters, we used 
niu = 0.364 GcV, which can be compared to the value of Mu(0) calculated here. We have 
also used either = 0.565 GeV [19-23] or = 0.585 GeV [18], values which may be 
compared to Ms(0). 

To proceed, we take A = 1.0 GeV, Gs = 13.30 GeV-^ k = 0.055 GeV^, m° = 0.0055 
GeV, m° = 0.130 GeV, // = 0.010 GeV and P = 20.0 GeV^ We then consider values of 
Gd ^ 0, Gd ^ -SOGs, Gd = -3OG5 and Go = -406*5. The results of our calculations 
arc given in Table 1. Recall that the function F{k) does not appear in our expression for 
the condensates. The calculation of the condensates includes the Pauli-Villars regulators, 
Cu{k^), Cdik"^) and C^(A;2), however. [See Eqs. (2.9)-(2.11).] In our calculation of the 
properties of the rj mesons [18] we had —Gd/Gs ~ 15 — 18, since we used Gd values in the 
range -180 GeV-^< Gg < -220 GeV"^ in that work. 

In order to specify a value oi Gd for this work, we note that a calculation based upon 
chiral perturbation theory yields {ss)/{uu) = 1.689 [24]. Inspection of Table I suggests that 
the values of Gd, other than Gd = 0, given in Table 1 arc acceptable. For Gd = —266 
GeV"^ we have M„(0) = 0.377 GeV and Ms{0) = 0.555 GeV, which are reasonably close to 
the phenomenological parameters = 0.364 GeV and nis — 0.565 GeV used in our earlier 
work [19-23]. 

It is worth noting that, in standard apphcation of the SU(3)-fiavor NJL model, one finds 
{ss)/{uu) ~ 1.1 [1], so that the results shown in Table I are encouraging, given that the 
value for {ss) / (uu) obtained using chiral perturbation theory is about 1.7 [24], as noted 
above. 

V. NUMERICAL RESULTS: 

MOMENTUM DEPENDENCE OF THE CONSTITUENT QUARK MASSES 

In Fig. 7 we show Mu{k), where k is the magnitude of the Euclidean momentum. The 
dashed line exhibits the result without the confining interaction (/t = 0). It is interesting 
to see that inclusion of confinement improves the shape of the curve when we compare our 
results to the lattice results shown in Figs. 1 and 2. We note that Mu{k) goes over to 
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Gd [ GeV-5] 


0.0 


-266.0 


-399.0 


-532.0 


M„(0) [GeV] 


0.334 


0.377 


0.396 


0.416 


M,(0) [GeV] 


0.538 


0.555 


0.564 


0.575 


{uu}^ [GeV] 


-0.207 


-0.215 


-0.217 


-0.220 


(ss)^ [GeV] 


-0.2605 


-0.261 


-0.261 


-0.261 


(ss) 


2.00 


1.80 


1.73 


1.68 


AuiQ) [GeV] 


0.447 


0.481 


0.496 


0.512 


S„(0) 


-0.335 


-0.276 


-0.253 


-0.233 


^,(0) [GeV] 


0.614 


0.628 


0.636 


0.645 


^3(0) 


-0.139 


-0.131 


-0.127 


-0.122 



TABLE I: Calculated values for the condensates and for A(0), 5(0), and M(0) are given for the up 
and strange quarks for four values of Gd- The parameters m° = 0.0055 GeV, = 0.130 GeV, 
K = 0.055 GeV2, (3 = 20 GeV«, /x = 0.010 GeV, A = 1.0 GeV, Gs = 13.30 GeV'^ were used. The 
values of k and /j, were fixed in earlier work [18-23]. Values of {uu) ~ {dd) ~ —(0.240 ±0.025 GeV)^ 
have been suggested [28], so we see that our calculated values are at, or near, the lower limit for 
that quantity. 

m° = 0.0055 GeV for large k. In Fig. 8 we show Bu{k). (Recall that Z„(A;) = [l-B^{k)]-\) 
We remark that Bu{k) — when k, — 0. In Figs. 9 and 10 we show Ms{k) and Bs{k), 
respectively. As expected, we find that Ms{k) goes over to m° = 0.130 GeV when k is large. 

VI. DENSITY DEPENDENCE OF THE QUARK MASS AND QUARK CONDEN- 
SATE 

As stated earlier, the behavior of the NJL model for finite values of the baryon density 
is an extensively explored topic [1, 25, 26], with particular recent emphasis on color super- 
conductivity [4-10]. In this section we explore the behavior of our model at finite baryon 
density. (It should be noted that the nature of the phase transition describing chiral sym- 
metry restoration at finite density is quite model dependent. For example, the inclusion of 
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FIG. 7: Values of M„(A;) are shown for the parameters m° = 0.0055 GeV, m° = 0.130 GeV, 
K = 0.055 GeV^, /x = 0.010 GeV, A = 1.0 GeV, /? = 20 GeV^ Gs = 13.3 GeV-^, and Go = -266 
GeV~^. The dashed hne shows the result without confinement (k = 0). 

current quark masses can change a strong first-order transition to a smooth second-order 
transition [26].) 

A comprehensive study of the thermodynamics of the three-flavor NJL model has been 
reported in Ref. [27]. There it is found that the up, down and strange quark masses are 
essentially constant up to the density where a first-order phase transition appears. At that 
point, the up and down quark masses drop from a value of about 380 MeV to about 30 MeV. 
That behavior differs from the behavior expected at low density. For example, we have the 
well-known relation between the value of the condensate and the baryon density of nuclear 
matter 

where ctjv is the pion-nucleon sigma term. This relation is valid to first-order in the density. 
It may be derived, in the case of nuclear matter, by writing 

{qq)p={qq)o + {N\qq\N)pB (6.2) 

and making use of the definition of the pion-nucleon sigma term, crjv, and the Gell-Mann- 
Oakes-Renner relation. If we put ctjv = 0.045 GeV, we have {qq)p/{qq)o — ^ — 0.273^^, 
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FIG. 8: Values of Bu{k) are shown. (See caption to Fig. 7.) 
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FIG. 9: Values of Ms{k) are shown. (See caption to Fig. 7.) 

where pb is in GeV^ units. For nuclear matter pB =(0.109 GeV)^, so we see that the 
condensate is reduced by about 35%, if we evaluate Eq. (6.1) at nuclear matter density. We 
can check whether the density dependence given by Eq. (6.1) is reproduced in our model, 
since it should not matter whether the scalar density of the background matter is generated 
by quarks in nucleons or by the presence of free quarks. In the former case, we may write. 
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FIG. 10: Values of Bs{k) are shown. (See caption to Fig. 7.) 
for the baryon density, 

- V (6.3) 

where the factor of 4 is arises from the product of the spin and isospin factors. In the case 
of quarks, we have 

"^ = ''"'{3) J Wf 

where, in this case, the factor of 4 again arises from the spin and isospin factor. (Both up 
and down quarks are present in equal numbers.) The color factor, N,, — 3, is cancelled by 
the baryon number of 1 /3 of each quark. 

We need to modify the equations for the quark self-energy to take into account the 
presence of the Fermi seas of up and down quarks whose Fermi momentum is kp- We take 
one Fermi sea to be composed of on- mass-shell up quarks with constituent mass Mu(0). The 
following term is then added to the equation for Au{k). 

4''W = -(2G.)iV.2/"g^/^(*-r) (6.5) 

1 

where Eu{k) = P + M2(0)] ' ' . The second factor of 2 in Eq. (6.5) reflects the spin degen- 
eracy. We note that M„(0) is density-dependent and could be written as Mu(0, p) in keeping 
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FIG. 11: The values of M{k,p) are shown for p/pNM = 0.26 [dotted Une], p/pNM = 0.52 [dashed 
Une], and p/pNM = 0.78 [dash-dot hne]. 

with the labeUing of Fig. 11, where Mu{k,p) is used. Also, if f{k — k') — 1, Au\k) would 
then represent —2GsPs, where ps is the scalar density associated with the up-quark Fermi 
sea. 

In Fig. 11 we show M{k, p) calculated for four values of p and in Fig. 12 wc show M(0) as 
a function of k\. Since the quarks in the Fermi sea are taken to be on- mass-shell, we would, 
in principle, require Mu{k,p) in Minkowski space. However, since kp — 0.268 GeV for the 
case of nuclear matter, only a very modest extrapolation of the curves shown in Fig. 11 is 
needed for the densities considered in this work. In Fig. 13 we show the value of the up quark 
condensate as a function of kp. (Note that kp = 19.2 x 10~^ GeV^ represents the density 
of nuclear matter.) It is seen that, for small values of the density, the density dependence 
of the condensate reproduces what is expected from Eq. (6.1). (If wc extrapolate the curve 
using a linear approximation, the condensate is reduced by about 30% at nuclear matter 
density.) 
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FIG. 12: Values of M(0) are shown as a function of kp. Note that ps = {2/3Tr'^)kp. 




FIG. 13: The values of the up quark condensate are given as a function of hp. Note that pB = 
(2/37r2)A:3. 

VII. DISCUSSION 

We have remarked earlier in this work that the large values of the condensate ratio 
(ss) I (uu) seen in Table I play a role in obtaining a good fit to the mixing angles of the 
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77(947) and 77/(958) mesons [18]. To understand this remark we note that the effective 
singlet-octet couphng constants for pseudoscalar states are [3] 

G^o = G'5-^(a + /3 + 7)^, (7.1) 

G^s = Gs- ^(7 -2a -2/5)^, (7.2) 



and 



Go'; = -^(27-«-/3)^, (7.3) 
where a — {uu), (3 — {dd) and 7 = (ss). We take a — (3, so that 

G£ = -^(7-")^. (7.4) 

If 7 = 1.7a, the result for is six times larger than when 7 = 1.1a. 

In addition to the effects of G^, singlet-octet mixing is induced by the quantity [18] 

Eos{k)^^[E^{k)-E,{k)], (7.5) 

where Ey_{k) — k'^ + , etc. It is found that, since and EQs{k) tend to cancel in our 
formalism, the significant singlet-octet mixing generated by Eg^ik) is reduced by the values 
of obtained for the larger value of the ratio (ss) / (uu) , with the result that we reproduce 
the values of the mixing angles found in other studies that make use of experimental data 
to obtain values for the mixing angles [18]. 

In our earlier work, which was carried out in Minkowski space, the values of m„ = rrid 
and nis were taken as parameters. Inspection of our figures which exhibit values of Mu{k) 
and Ms{k) suggests that an extrapolation into Minkowski space may be made if k"^ is not 
too large. The fact that Mu(0) and Ms(0) are close to our phenomenological parameters for 
Gd = —266 GeV~^ is encouraging and suggests that some support for our choice of quark 
mass parameters may be found in our Euclidean-space analysis. 

The full consequences of separating the specification of the nonlocality of the quark inter- 
action from the choice of the regulator of the theory should be explored more fully. Although 
that feature of our model introduces greater flexibility, that comes with the disadvantage of 
having to introduce other parameters in the model. We have made only limited variation of 
the form of the nonlocality and the regulator. For further applications it may be of interest 
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to explore a more comprehensive parameter variation. It is also necessary to extend the 
calculations reported in Figs. 11-13 to larger values of the density than those considered 
here. That step will require more complex methods for solving our nonlinear equations for 
the self-energy than the simple iteration scheme we have used thus far. 

Our work may be compared to that of Alkofer, Watson and Weigel [29] who have solved 
the Schwinger-Dyson equation using a gluon propagator whose low-momentum behavior 
is enhanced by a Gaussian function. (That modification requires the introduction of two 
phenomenological parameters [30].) The behavior found for A{k) and B{k) in Euclidean 
space is similar to that obtained in this work. (See Fig. 1 of Ref. [29].) Those authors 
also solve the Bethe-Salpeter equation to obtain the properties of various qq mesons with 
generally satisfactory results. It is of interest to note that the Minkowski space solution for 
A{k) and B{k) is such that the quark can go on- mass-shell. That feature may be related to 
our work [18-23] in which we use on-mass-shell quarks with masses = = 0.364 GeV 
and TUs = 0.565 GeV (or 0.585 GeV [18]) when solving the Bethe-Salpeter equation in our 
study of qq mesons. 
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